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ABSTRACT 

If X-ray flashes (XRFs) and X-ray rich Gamma-ray Bursts(XRRGs) have the same origin 
with Gamma-ray Bursts (GRBs) but are viewed from larger angles of structured jets, their 
early afterglows may differ from those of GRBs. When the ultra-relativistic outflow interact 
with the surrounding medium, there are two shocks formed, one is a forward shock, the 
other is a reverse shock. In this paper we calculate numerically the early afterglow powered 
by uniform jet, Gaussian jet and power-law jet in the forward-reverse shock scenario. A 
set of differential equations are used to govern the dynamical evolution and synchrotron 
self-Gompton effect has been taken into account to calculate the emission. In uniform jets, 
the very early afterglows of XRRGs and XRFs are significantly lower than GRBs and the 
observed peak times of RS emission are longer in interstellar medium environment. The 
RS components in XRRGs and XRFs are difficult to be detected. But in stellar wind, the 
reduce of very early flux and the delay of RS peak time are not so remarkable. In nonuniform 
jet(Gaussian jet and power-law jet), where there are emission materials on the line of sight, 
the very early light curve resembles isotropic-equivalent ejecta in general although the RS 
flux decay index shows notable deviation if the RS is relativistic(in stellar wind). 
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1 INTRODUCTION 

Late time Gamma-ray Burst (GRB) afterglow emission has 
been well observed and studied since the first detection in 
1997 ivan Paradiis et al ] ll997l: iGosta et al In contrast, 

afterglow shortly after or during the main GRB is still quite 
uncertain because of the lack of observation. This early af¬ 
terglow may provide important information about initial pa¬ 
rameters of the burst and shed light on explosion mechanism. 
SWIFT is able to observe multi-wavelength radiation rapidly 
after 7 ray triggering. With new observations, the early stage 
should be investigated more profoundly. In the standard bary- 
onic GRB fireball model, prompt GRB emission is produced 
by collisions of internal shocks while afterglow comes from in¬ 
teractions between burst-ejected materials and circum-burst 
medium. After the internal shock phase, as the fireball is decel¬ 
erated b^_th££ircunibmstjnedlurnj_usunlRa_£mr of shocks de¬ 
velop iMfezaros fc R,eejll997l : ISari fc Piranlll999ll . One is a for¬ 
ward shock (FS) expands outward to heat the external medium, 
the other is a reverse shock (RS) propagates into the ejecta 
and heats this cold initial shell. The central energy is gradu¬ 
ally transmitted to outer medium through shocked mediums 
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and shocked ejecta. The early afterglow studied in this paper is 
during this transition stage. 

iMeszaros fc Ree'^ (Il99i1) pointed out RS should emit 
detectable optical synchrotron photons. This prediction 
has been confirmed by the optical flash and t h e ra- 
dio flare detected in GRB 990123^Akerlo^^lJ_|ffi9^ 

g a mLeU^ [l999|; JSari fc PiranI Il99^ | Mfezd^s fc R.eed| 
~~ l^[ki^fcpiraM I20 o3l ~ T he later more detailed in- 
vestiga tion suggests that the RS region may be magne- 
tized _l|Faj 2 _et_aJj 200 j|_ JZhang. Kobavashi fc Meszaro^ 1200,31 : 
|Panaitesci^^Gnnaill2004l) . The RS emi ssion of GRBs in stel¬ 
lar wind medium has been discusse d llChevalier fc Lil 120001: 
1 ^ et al ] 1200.31: IZou. Wu fc Daill200!^^ . The RS emission pow¬ 
ered by the magnetized outflow o r neutron-fed outflow has 
been investigated in some detail _l|Faji ^_^i_fcAVan3 l2004bt 
IZhang fc Kobavashil200,‘Tl : iFan. Zhang fc Weill200^al) . The pair- 
rich re verse shock emission has been discus sed by iLi et al.l 
i2003ari and lMcMahon. Kumar fc PiranI ||200^. Unfortunately, 
up to now, there are only a few additional candidates. They 


are GRB 021211 JFox ] l200,3l: iLi et aljEio^ IWc. 
Kum^^^^anaite^^^^T^ GRB 041219a (Blak^^^I 


^^^han T^^^l^OOfiT and GRB 050525a 


2003; 


2005; 


200 ," : 


Shao fc Daill200!Ttl . Whether the lack of RS emission events is 


due to observational limits or theoretical problems may be set¬ 
tled by further observation. 
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X-ray Flash (XRF) is an interesting phenomenon 
which resembles GRB. It is similar to GRB in many 
aspects except its lower peak energy (~10keV) and 
lower is otropic energy 10^^ergs ). Sever al models are 

g ose djjDerme^^hlan^&Bottchea [l999|; _|Heiseet^^ 

J20^ iBaraaud et al\ l200S ^ 
3. Donaghv fc Grazianil l200?^^ among which the struc¬ 
tured jet model is widely discussed. When a jet is ejected 
and propagating, observers on different viewing angles will 
see different phenomena. GRBs are observed near the jet 
center. XRFs may be detected at the edge of the jet. Be¬ 
tween them are the X-ray rich GRBs (XRRGs) whose peak 
energies and isotropic energies are between those of GRBs 
and XRFs. Possible jet models for the X RFs and XRRGs 
inclu d e the off-beam uniform jet mode l llloka fc Nakamuril 
1200ll : IWrniazak i, loka fc Nakarnural IjPOd), the Gaussian jet 
model l ^ang fc MeszaroslbooJTLlovXRonnin g. Dai fc Zhand 

I 2 OO 4 I : the power-law 

jet model iMfezdro^^n ~ ^9^ I.Tin fc Weil l2004l . The very 
ear ly afterglow powered by st ructured jets has been calculated 
by iFan. Wei fc Wan3 (l2004a^ analytically. In this work, we 
present our numerical results. 

The dynamical evolution is described in Section 2. The ra¬ 
diation calculation procedure is shown in Section 3. Afterglows 
from jets are prsented in Section 4. In section 5, our results are 
summarized with some discussions. 


2 DYNAMICAL EVOLUTION 

Consider a homogenous cold shell ejected from central engine 
with isotropic energy Eiao and width Aq. It expands relativis- 
tically with Lorenz factor r) and contains mass rriej = Eisolr\<?. 
When a FS and a RS are formed, there is a shocked region con¬ 
sists of tow parts: the shocked medium and the shocked shell 
separated by the contact discontinuity surface. The shocked re¬ 
gion is assumed to have homogenous bulk velocity and energy 
density, that is, 72 = 73, 62 = 63 (hereafter the subscript 2 
denotes the shocked medium, 3 denotes the shocked shell and 
4 denotes the cold ejecta). However, the number density and 
random particle velocity in the comoving frame are discontinu¬ 
ous at the contact discontinuous surface. The random Lorentz 
factor of the shocked protons are (72 — 1 ) ~ 72 and (734 — 1 ) in 
regions 2 and 3 , respectively (iBlandford fc Mckeelll 97 ^ . where 
734 is the relative Lorentz factor of the shocked shell to the cold 
shell. 

Kinetic energies of the shocked medium, the shocked shell 
and the unshocked shell are respectively: 

El = (72 — l)m2C^ + (1 — £2)72^2, 

E3 = (73 — l)m 3 C^ + (1 — £3)73^3, 

E4 = {rj - l){mej - mz)c, 

where U2 = (72 —l)m2C^ and U3 = (734 —l)m-3C^ are internal en¬ 
ergies in the comoving frame. If a fraction e of the fresh shocked 
material energy is radiated, the radiated thermal energies of the 
shocked medium and shocked shell are £272(72 — l)dm2C^ and 
£373(734 — Vjdmsc?■ With energy conservation, we have 


If GRBs are located in interstellar medium(ISM), the 
number density of the external medium ni is constant. 
If GRBs are born in stellar wind, ni = AR~^, where 
A = M/A-KmpV-u, — 3 X 10®®and A, = 

(M/1 X 10~®M^t/r~^ln,„ /10®fcm/si~^ iChevalier fc Lil 1200(1) . 
The swept-up mass of FS evolves as 

—— = 4nR mmp. (2) 

dH. 

The comoving number density of the ejecta 714 = 
Mej/[ 47 rmpR^A? 7 ], where A = max{Ao, R/r;^) is the width of 
shell considering spreading effect. When R is smaller than the 
spreading radius Rs = Aq??^, A ~ Aq. If R > Rs, A ~ Rlv^- 
The swept-up mass of RS evolves as 


= 4tvR'^{P 4 - PRs)rjn4mp, (3) 

where Prs ~ ’ which is the Lorentz factor of RS 

jFan. Wei fc War^l2()04bli . From equations (EP, it is found 
that 72 evolves as 


d'y2 

dR 


= -477 


Q 

R’ 


where 


(4) 


Q = {'ll - l)nimp -I- (72734 - r]){r]n 4 mp){P 4 - Prs), 

P = rn2 + m3 + {1 — £2) (272 — l)m2 + (1 — £3) (734 — l)m3 

+ (1 - £3)72m3[77(l - P 2 P 4 ) - 

72/22 

The overall dynamical evolution of RS-FS can be obtained by 
solving equations lIM combined with 


dR = 


P 2 cdt 
1 — P 2 1 + 2 


(5) 


where t is the time measured by the observer and 2 is the red- 
shift. 

Define rjs to be the electron radiative efficiency, that is, a 
fraction of the electron energy is radiated. In the fast cooling 
stage, electrons cool imm ediately and rie = 1 . In the slow cooling 
stage, 77e = (7 m /ycY ^ llSari fc EsirJl200ili . where 7m and 7 c 
are the minimum and the cooling Lorentz factor of the shocked 
electrons, respectively (see detail in the next section). Usually it 
is assumed that electrons share a fraction Ce of internal energy. 
Then the radiative efficiency is 

£ = £e?7e = £e min[l, (7m/7c)'’~^]. ( 6 ) 


e = 0 corresponds to the adiabatic case while e = 1 corresponds 
to t he fully radiat i ve ca se. 

ISari fc PirarJ lll99?jl introduced ^ to de¬ 

scribe the strength of the RS, where I = [ 3 Riso/( 47 rnimpC®)]'^/® 
is the Sedov length. If ^ < 1 at the RS crossing time (tx), the 
RS is relativistic (RRS). If ^ > 1 at fx, the RS is non-relativistic 
(NRS). Our results are consist with them. The scaling laws are 
exactly the same with analytical results though differences exist 
in absolute values. In our numerical results, it is found out that 
734 — 1 can be well fitted by 


lg(734 - 1) = -0.64991 - 1.64406 IgC 

-0.59494(lg - 0.13656(lg 0^. (7) 


d{E 2 + E 3 + E 4 ) = —£ 272(72 — l)dm 2 C 

-£373(734 - l)dm3C^. (1) 


for 0.01 < C < 10 . 

After the RS crosses the shell, the shocked medium and the 
shocked ejecta evolve somewhat independently. The FS con¬ 
tinues to spread outwards, collect external medium and the 
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shocked medium is gradually shaped into BM self-similar pro- 
hle. The energy conservation still applies: 

di ?2 =-£ 272(72 - (t>tx). ( 8 ) 

Then the forward shock evolution can be solved together with 
equation(|2l, 0 and 0. At the same time, the shocked ejecta 
ceases to increase, he. , m 3 = msj = const for t > tx- 
iKobovashi fc Sail Il2n0fil showed that its evolution consists ap¬ 
proximately with BM solution, that is, 73 oc {k = 0 

applies to ISM and A: = 2 to wind). They also inferred that the 
shocked ejecta spreads with the speed of light in the comoving 
frame and evolves adiabatically. 


3 RADIATION 


Both FS and RS heat the cold materials to higher tempera¬ 
tures, accelerate protons and electrons as well as generate ran¬ 
dom magnetic field. Here we calculate the synchrotron radiation 
and the synchrotron self-Compton cooling of the shocked elec¬ 
trons. The internal energy density of FS shocked medium 62 = 
do-oni m„c^ (lBlandford fc Mckeelll97(il . For the ejecta shocked 

_ 4(3+g) 

by RS, 63 = 62 before tx , otherwise 63 ext 3 ( 1 + 29 ) (^g — r j2 — k). 
Ce and eb are defined to be the faction of internal energy oc¬ 
cupied by electrons and magnetic held. It is easy to get that 
magnetic held B = VSneBC. As to electrons, in the absence 
of energy loss, they are assumed to be shocked to a power-law 
distribution dNs/d'y^ oc 7T^(7m < 7e < 7 m) with minimum 
Lorentz factor 


'ym ,2 


£e(72 


^ mp(p- 2) 

mjp - 1 ) 


+ 1? 


'TttIjS 


£e(734 - 


mp{p- 2 ) 

me{p- 1) 


+ Ij 


and m aximum Lorentz factor jm ~ 10®(B/1G)’ ^^^ dPai fc LrJ 
I2OOIII . He re e and B are all measured in the coming frame. 
Following ISari et ^ il998li , we introduce the cooling Lorentz 
factor 7 c = to describe the synchrotron radiation 

loss of the shocked electrons. Note that 7 c ,3 = 7 c ,2 before tx- 
The act ual electron distrib u tion sh ould be given as following 
(see also iFan. Wei fc Wand (l2004lJH : (i) for 7 c ^ 7 m, i.e., the 
fast cooling case 


dNe _ n i ^ 7e ^ 7iii), 

d7e ^ 1 , ( 7 ™ < 7e < 7 m), 


(9) 


where 

Gi = [— - — 

7c P 7m 


■p —1 —P 

7m 7 m 1-1 


Ntc 


where Atot is the total number of radiating electrons involved; 
(ii) for 7m < 7c ^ 7 m, i.e., the slow cooling case 



7c 7 . 


-(p+i) 


where 


-l-P ryi-p 

Im I c 

p-i p(p-i) 


(7m ^ 7e 7c), 
(7c < 7e ^ 7 m), 


7 c 7 m^ 

p 


r^Atot. 


( 10 ) 


In the comoving frame, the synchrotron radiation of the 
two-segment distributed electrons can be well described by 
power-law spectrum consisted of se veral segments. The sp ec- 
trum peaks at min(i 2 m, 1 /') with hux dWiiers fc Galamal[T999^ 


F' = 

m.n.er. ^ '1 


VSe^BNtot 


rrieC^ 


(11) 


where "Fp is a function of p(for p ~ 2.5,4>p ~ 0.60). The 
breaking frequencies corresponding to 7 m and 7 c are vln = 
Sjf^eBf (A-Kniec) and p = 37 ceR/( 47 rmec). The absorption fre¬ 
quency n'a , under which the synchrotr on self-absorption can not 
be ignored, is calculated accordingly IWu et fflZl feOO.'ii . These 
three breaking frequencies divide the spectrum into four seg¬ 
ments. The spectrum indices are, from low frequency to high fre¬ 
quency, (i)for 1 /', [2, 1/3, (l-p)/2, -p/2]; (ii)for i// < 

12 / < [2,1/3, -1/2,-p/2]; (iii)for v'^<v',< P, [2,5/2, (1- 

P)/2,-p/2]; (iv)for p < < vV [2, 5/2,-1/2,-p/2]; (v)for 

v'm<v'o< Pa and < 1 //, [2, 5/2, 5/2, -p/2]. 

We assume that the synchrotron power is radiated 
isotropically in the comoving frame, ^ 

T he angular distribution of power in the ob server’s frame 
is te.vbicki fc LightmarJll976l : lHuang et ^ 11 ^ 


1 -b z dF'(iy') 


dF(iy) _ 
dfl ^ 73 ( 1 -/3p)3 dfl' 

where 


1 + z F'(F) 

73(1 —/?p)3 47r 


( 12 ) 


(l-h z)j(l- p/d)' 

Then the observed flux density at frequency iz is 


S. - j 


dF(n) A 


dn 




( 13 ) 


(14) 


where A is the area of our detector and Bl is the luminosity 
distance (we assume Hq = 65km s“^ Mpc“^, Qm = 0.3, Ha = 
0.7). 

After the RS crosses the shell, no new electrons are accel¬ 
erated in the shocked shell anymore and all the electrons cool 
at the same rate caused by the adiabatic expansion of the fluid. 
No electrons exist above 7 c, so the flux above zz/ drops expo¬ 
nentially. 

At a particular time t, the photons received by the ob- 
server comes from a n equal arriving time surface determined by 
iHuang et al 11 ^ 


t — (1 -|- z) 


1-/3 cos 0 
/dc 


dR = const.. 


(15) 


where 0 is the angle between the radiation region and the view¬ 
ing line in the burst frame. 

We also consider synchrotron self-Compton (SSC) effect. It 
will reduce the cooling Lorentz factor to 


7c = 


Ic 

1 + Y’ 


(16) 


where Y is the Compton par ameter, expressed as jSari fc EsirJ 
l200ll:IWang. Dai fc L)Jl200lll 


Y _ Bias _ PetalCB _ —1 + \/l + dpeteJtB , . 

^ Leya ^ l-bT “ 2 ■ 

Take it into eauation I16t , 7 c can be solved numerically. At early 
time, the forward shocked material is in the fast cooling phase 
and the cooling frequency will be reduced signihcantly due to 
the SSC radiation. As to the RS, whether it is in the fast or slow 
cooling phase is parameters dependent. However, even if it is in 
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slow cooling phase, the cooling frequency will be suppressed by 
SSC if te/^B S> 1, which is the ordinary case. 


4 NUMERICAL RESULTS 

In this section the early afterglow from isotropic ejecta is pre¬ 
sented at first. Then we calculate afterglows powered by differ¬ 
ent jets. Three kinds of jets are considered. They are the uni¬ 
form jet, the Gaussian jet, and the power-law jet. The sideways 
expansion of these jets is ignored in our calculation. 

The obstacle we encountered in the current work is the 
poorly known initial Lorentz factor of these jets, especially for 
the structured ones (The Gaussian jet and power-law jet). By 
now, the spectrum of most XRFs are non-thermal. The emitting 
region should be optical thin, i.e., the optical thickness r = 
rieaTR' < 1, where R' is the comoving width of the emitting 
region, comoving number density 

of electrons contained in the emitting region. L is the isotropic 
luminosity of the burst, St is the typical variability timescale 
of the burst lightcurve. Assuming that XRFs are powered by 
internal shocks, t he optical thin (r < 1) condition yields (e.g., 
iFan fc Weil (1200, 

rj ^ 5^2^®. (18) 

iKumar fc Grano^ i200,^l have performed a numerical inves¬ 
tigation on the hydrodynamical evolution of a Gaussian jet by 
assuming ri{9) is also Gaussian. In the Monte Carlo simulation 
of Zhang et al. (2004), ri{9) has been taken as a free function 
but to find that a fiat or only small fluctuation angular distribu¬ 
tion of Lorentz factor is indeed required. Therefore, and partly 
for convenience, we take rj = 300 for off-axis uniform jet, the 
Gaussian jet, and the power-law jet. 

For each jet model, light curves with typical parameters 
from three viewing angles are presented, which may correspond 
to GRBs, XRRGs and XRFs respectively. The light curves differ 
from each other although there are some universal properties. 
Viewing from the center of a jet, the light curve is approximately 
the same as that powered by an isotropic ejecta at the early time 
and a break occurs when the emission from the edge comes into 
view. When we observe off-axis, the flux is lower at the early 
time but gradually converge into the light curve from the center 
when the whole jet can be seen. 

4.1 Isotropic ejecta 

The upper panel of Fig.l shows the early R band light curve 
from isotropic ejecta in ISM case with typical parameters. The 
RS dominates the very early R band light curve, peaks at about 
14mag. The RS is non-relativistic at tx with < 

vr < . Therefore tx ~ 45s is greater than the prompt du¬ 

ration Tdur = (1 + -z)-^ ~ 20s. The R band arising index is 
^2.5, which is shallower than the analytical estimation of NRS 
iKobovashl200ntlFan etolj gOO^lMt^ v aries w ith time as well as 
with initial nara.meters jNak^^^ Piranll2nn4l . After tx , the de¬ 
cay index is ~ 2 until = ur. Then a break occurs followed 
by faster decay. The flux after the break comes from equal ar¬ 
riving time surface emitted at lower radius. Since Vc is reduced 
by SSG, it reaches the R band and the RS flux encounters faster 
decay shortly after tx - Due to the rapid decay of RS emission, 
the later emission is mainly contributed by FS. It is interest¬ 
ing to note that the late FS emission bump vanishes when the 



10 ° 10 ’ 10 ° 10 ° 10 ’ 10 ° 10 ° 
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Figure 1. R band afterglow from isotropic ejecta with typical pa¬ 
rameters. Upper panel is the ISM case and lower is the wind case. 
Thin lines are the RS emission and thick lines are the total emission. 
SSC effect is ignored in the dashed lines and taken into account in 
the solid lines. In the ISM case, the shape of the light curve changed 
since , 1 /^^ ~ ub around tx - The RS shows a more rapid decay 
and the bump caused by the FS vanishes. In the wind case, both 
RS and FS flux are reduced since max(i/,?‘®, < Rr- Parameters: 

rj = 300, Aq = 3 X 10^'^cm, Eiso = 1 x 10®®ergs, ni = lcm“® for 
ISM and r) = 300, Aq = 1 x Eiso = 5 x 10®®ergs, A* = 1.0 

for wind, Ce = 0.3, es = 0.01, p = 2.5, z = 1- 


SSG effect has been taken into account. The physical reason is 
that with SSC effect, the FS flux peak time appears earlier at 
rr = during fast cooling. Since r^^ oc oc t~^^^ 

and ~ const, the decay index following the peak flux is 

quite flat(-l/4) and change to —(3p — 2)/4 after rr — R^^- 
In stellar wind environment, the RS also dominates the 
very early R band light curve with typical parameters, which is 
shown in the lower panel of Fig.l. The peak flux is about 13mag. 
The RS is relativistic at tx with r^^ < r^^ < Rr < Rm^- 
tx ~ 40s coincides with Tdur ~ 65s. The R band arising index 
is ~ 0.5, which is much flatter than that in the ISM case. After 
tx, the flux decay with -2.5, which is contributed to emission 
from lower radius on equal arrival time surface. Since both RS 
and FS are in the fast cooling case and max{i/,f‘®, < rr, 

the R band flux is reduced by a factor 1 + Y ~ l-|--\/ee/es(about 
2 magnitude taking = 30) when the SSC effect has been 

taken into account. 
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4.2 Uniform jet 

Uniform jet is a simple model which assumes that all the ejected 
materials are confined to a collimated uniform cone. GRBs are 
detected within the cone (on-beam). XRFs and XRRGs are 
viewed out of the jet (off beam). If our viewing angle ©t, is 
slightly larger than the jet opening angle 9jst, the observed fre¬ 
quency 

^off — CLiyonj 

where a « [1 -|- 7^(A0)^]“^ (7 is the bulk Lorentz factor of 
the jet) and A© = 0„ — 6 jet- Taking the peak energy of XRFs 
Ep.xRF ~ O.lEp^GRB, XRFs should be observed at AQxrf ~ 
3/7. The observed flux can be es timated by an empirical formula 
(e.g.. lFan. Wei fc Wand ijioM^)) 

3 

Evoffi^^Jobs) ~ -^E,^^„{0,t), ( 20 ) 

where dtobs = dt/a. It implies that if 7 > 1/A©,the observed 
flux is dramatically lower and the observed time is longer. The 
early afterglow powered by the uniform jets has been presented 
in Fig.|^ 

In the ISM case, the flux increases rapidly at the early time. 
The RS component reaches its peak flux at a time longer than 
viewing on-beam. The peak flux of RS emission for A© = 3/(277) 
(i.e., the XRRG) is about 17.5mag at about 100s. It is bright 
enough to be detected. The rising index of RS is about 3.5. The 
peak flux of the RS emission for A© = 3/7 (i.e., the XRF) 
is about 21mag at about 250s. It is hard to be detected by 
UVOT on SWIFT or ground telescope at such early time. The 
rising index is about 4.2. After tx, the observed FS flux keeps 
rising for some while, peaks several hours later. The peak FS 
emission flux of XRRG is about 17mag and 17.5mag for XRF. 
Therefore in XRFs, the RS component can hardly be seen and 
the later FS component around Ihour is much more luminous. 
In XRRG, the RS emission is detectable but the RS component 
is not significant. The light curve is nearly flat for t <1 hour. 

In the wind case, the off-beam RS very early flux increases 
at the same rate as in the on-beam case since both 7 and a 
are constant, but the flux is lower. The peak time of off-beam 
RS is slightly longer, ~ 50s in XRF. The XRF peak flux is 
about 15mag and significantly higher than the FS. After tx, as 
the Lorentz factor decreases, the light curve gradually converge 
into the light curve from the center. 

4.3 Nonuniform jet 

Different from uniform jets, there are emission materials on 
the line of sight in nonuniform jets(Gaussian jets and power- 
law jets). XRFs and XRRGs are viewed within the jets but 
have lower isotropic-equivalent energy than GRBs. The ob¬ 
served early emission is mainly contributed from materials on 
the line of sight. However, we find that the light curve show 
notable deviation if the RS is relativistic. 

4-3.1 Gaussian jet 

The jet energy in a Gaussian jet distributed with angular as 

E{e) = 0 s; 0 (21) 

With typical value that Eiao.xRF ~ the viewing 

angle of XRF should be Qv,xrf ~ 3.04^o. We also calculate 



Figure 2. R band afterglow from uniform jet. Upper panel is the ISM 
case and lower is the wind case. Thin lines are the RS emission and 
thick lines are the total emission. A© = 0, 3/27;, 3/?7 in solid lines, 
dashed lines and dotted lines respectively. Parameters: 9jet = 0.1, 
other parameters are the same as in Fig.l. 


the light curve with Ov,xrrg ~ 2.15^o which may correspond 
to the XRRG(taking Eiso,xRRG = 10 ^Eiso ,grb)- The R band 
afterglow from Gaussian jet is shown in Fig.3. 

In the ISM case, the off-axis RS flux is exactly the same 
with isotropic emission taking Eiso = E{Ov)- The RS is non- 
relativistic and the Lorentz factor is high at early time. So the 
jet structure can not change the light curve significantly. As ©„ 
increases, the light curve has (i) shorter (ii)lower peak RS 
emission flux, about 17mag for XRRG and 19.5mag for XRF; 
(iii) flatter rising before tx', (iv) longer ti{c(the time ur = Uc), 
therefore longer time interval for the ~ —2 delay after tx- 

In stellar wind, the off-axis RS flux is slightly different 
with the isotropic emission taking Eiso = i5(©^). The RS is 
relativistic and the Lorentz factor is sufficiently low to reflect 
the jet structure. The isotropic emission with decreasing Eiso 
from GRB to XRF has (i)nearly constant tx, constant rising 
index 0.5 before tx and decay index -2.5 after tx ; (ii)lower peak 
flux, about 16.5mag for XRRG and 19mag for XRF. What ac¬ 
tual light curve differs from isotropic one is the decay index of 
RS flux. Viewing from the center, RS flux decays faster. With 
increasing viewing angle, the start time of the faster decay in¬ 
creases, and from t x to this time, the flux is flatter than isotropic 
emission. The surplus emission comes from the center of the jet 
with 6 < Oo- 
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Figure 3. R band afterglow from Gaussian jet. Upper panel is the 
ISM case and lower is the wind case. Thin lines are the RS emission 
and thick lines are the total emission, ©u = 0, 2.15^0? 3.04^0 solid 
lines, dashed lines and dotted lines respectively. Parameters: Eq = 
1X for ISM and Eq = 5 x 10^^ergs for wind, Oq = 0.05,0jet = 

0.2, other parameters are the same as in Fig.l. 


Figure 4. R band afterglow from power-law jet. Upper panel is the 
ISM case and lower is the wind case. Thin lines are the RS emission 
and thick lines are the total emission. ©^; = 0, 3.160c, 100c in solid 
lines, dashed lines and dotted lines respectively. Parameters: Ec = 
1 X 10^^ergs for ISM and Ec = 5 x lO^^ergs for wind, 0c = 0.015, 
Ojet = 0.2, other parameters are the same as in Fig.l. 


4 . 3.2 Power-law jet 


In this model, the jet energy distribution as a function of angular 
is as follows 


E{6) 


Ec 0 ^ 0 < 0c 

Ec{6l6c)~^ 6c ^ 9 ^ Ojet- 


( 22 ) 


EiBo,XRF ~ 10 ^Eiso,GRB and Eiao,XRRG = 10 ^Eiao,GRB Cor¬ 
respond to Qv,xrf ~ 100c and Qv,xrrg ~ 3.160c respectively. 
Early afterglow from power-law jet is shown in Fig.4. 

In the ISM case, the RS is also exactly the same with 
isotropic emission taking Eiao = i5(0„), as found in the Gaus¬ 
sian jet model. In the wind case, RS flux also resembles that in 
the Gaussian jet except the decay index departure from isotropic 
emission is more prominent. 


5 CONCLUSION AND DISCUSSION 

So far, there are several RS emission candidates reported in 
SWIFT era. One is GRB 041219a detected by INTEGRAL. 
iFan. Zhang fc i2005hr) fitted its early IR afterglow with 
RS-FS emission model and found the RS region is magnetized. 
GRB 050525a may also be an possible candidate. IShno fc Pail 
^2006^ used RS-FS in standard scenario to fit the early optical 
bump and the parameters are reasonable. 


However, among all the bursts targeted in optical band 
within few minutes after the prompt emission, only a small frac¬ 
tion of them have likely RS emission. This is conflicted with the 
theoretical estimation. It may be caused by the extinctions of 
the host galaxy. And as demonstrated in this paper, the over¬ 
estimation of the Vm, Re may lead to a dimmer RS radiation. 
Up to now, the very early optical afterglows of XRFs have not 
been well detected, future observation may help us to modify 
the present theory. 

In this paper, we calculate the early afterglow powered by 
various kinds of jets numerically. Dynamical evolution is solved 
from a set of differen tial equations. This i s differ ent from the an¬ 
alytical treatment of iFan. Wei fc Wand (l2004a^ and gives simi¬ 
lar but more exactly results. We find that the most unprecise es¬ 
timation comes from 734 — 1, which results in the overestimation 
of Vrn and too rapid increasement of Vrn before tx in previous 
treatment. This is especially significant in the ISM case since 
the RS is mildly relativistic in typical parameters taken here. At 
the same time, Vc was overestimated previously because of the 
ignoring of SSC effect. Considering these two factors, the peak 
flux of reverse shock should be dimmer, which may help us to 
explain the lack of detection of optical flashes in most GRBs. It 
also needs to be pointed out that electrons may be cooled not 
only by SSC process but also by external photons) e.g., prompt 
emission) through inverse Compton scattering. Then Vc may 
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be reduced further and lead to even fainter very early optical 
radiation. 

Early afterglow from jets, both uniform and structured, are 
calculated. We find that the early afterglow varies signihcantly 
with different viewing angles and is dependent on the jet struc¬ 
ture. 

SWIFT XRT detection hnd that the early X-ray flare 
may be an comm on characteristic of GRB X-ray afterglow. 
iFan fc Weil ^200!^^ suggested that the RS synchrotron emission 
hardly can produce the very early X-ray flares. Our calcula¬ 
tion confirms their results. We also tried to simulate it with RS 
SSC emission within a large parameters space but failed. The 
RS SSC X-ray emission is always lower and generally far lower 
than FS synchrotron emission. 
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